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Abstract 

We investigate classical and quantum correlations of a quantum field in the inflationary universe 
using a particle detector model. By considering the entanglement and correlations between two 
I comoving detectors interacting with a scalar field, we find that the entanglement between the 

CN ' detectors becomes zero after their physical separation exceeds the Hubble horizon. Furthermore, 

W)' the quantum discord, which is defined as the quantum part of total correlation, approaches zero 

on the super horizon scale. These behaviors support the appearance of the classical nature of the 
quantum fluctuation generated during the inflationary era. 
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I. INTRODUCTION 



According to the inflationary scenario of cosmology, all structure in the Universe can be 
traced back to primordial quantum fluctuations during an accelerated expanding phase of 
the very early universe. Short wavelength quantum fluctuations generated during inflation 
are considered to lose quantum nature when their wavelengths exceed the Hubble horizon 
length. Then, the statistical property of generated fluctuations can be represented by clas- 
sical distribution functions. This is the assumption of the quantum to classical transition of 
quantum fluctuations generated by the inflation. As the structure of the present Universe is 
classical objects, we must explain or understand how this transition occurred and how the 
quantum fluctuations changed to classical fluctuations. 

We have investigated this problem from the viewpoint of entanglement 0]. For quan- 
tum fluctuations behaving in the classical way, quantum expectation values of any operators 
must be calculable using appropriate classical distribution functions. Restricting to two 
point correlation functions of quantum operators, this condition is equivalent to the separa- 
bility of a bipartite state. For a separable state, the entanglement is zero and there exists 
a positive normalizable P function [3] and this function can play the role of classical distri- 
bution function. In our previous study, we defined two spatially separated regions in the 
inflationary universe and investigated the bipartite entanglement between these regions. We 
found that the entanglement between these two regions becomes zero after their physical 
separation exceeds the Hubble horizon. This behavior of the bipartite entanglement confirms 
our expectation that the long wavelength quantum fluctuations during inflation behave as 
classical fluctuations and can become seed fluctuations for the structure formation in the 
Universe. Our previous analysis concerning the entanglement of quantum fluctuations in the 
inflationary universe relies on the separability criterion for continuous bipartite systems 0, Isf 
of which dynamical variables are continuous. The applicability of this criterion is limited 
to systems with Gaussian states: the wave function or the density matrix of the system 
is represented in a form of Gaussian function. Thus, we cannot say anything about the 
entanglement for the system with non-Gaussian states such as excited states and thermal 
states. Furthermore, from a viewpoint of observation or measurement, information on quan- 
tum fluctuations can be extracted via interaction between quantum fields and measurement 
devices. Hence, it is more natural to consider a setup in which the entanglement of quantum 
fields is probed using detectors. 

Following this direction, we consider particle detectors [gI, 0] with two internal energy 
levels interacting with a scalar field in this paper. By preparing two spatially separated 
equivalent detectors interacting with the scalar field, we can extract the information on 
entanglement of the scalar field by evaluating the entanglement between these two detectors. 
As a pair of such detectors is a two-qubit system, we have the necessity and sufficient 
condition for entanglement of this system js], [9|. Using this setup, B. Reznik el al. [13, 



111 ] studied the entanglement of the Minkowski vacuum. They showed that an initially 
nonentangled pair of detectors evolved to an entangled state through interaction with the 
scalar field. As the entanglement cannot be created by local operations, this implies that 
the entanglement of the quantum field is transferred to a pair of detectors. M. Cliche 
and A. Kempf [l2| constructed the information-theoretic quantum channel using this setup 
and evaluated the classical and quantum channel capacities as a function of the spacetime 



separation. G. V. Steeg and N. C. Menicucci [13| investigated the entanglement between 



detectors in de Sitter spacetime and they concluded that the conformal vacuum state of the 
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massless scalar field can be discriminated from the thermal state using the measurement of 
entanglement. 

In this paper, we investigate the entanglement structure of the quantum field in the 
expanding universe using the particle detector model. We also consider correlations between 
detectors and explore the relation between classical and quantum parts of correlations. This 
paper is organized as follows. In Sec. II, we present our setup of a detector system. Then, in 
Sec. Ill, we review entanglement measure (negativity) and classical and quantum correlations 
for a two-qubit system. In Sec. IV, we calculate entanglement and correlations for quantum 
fields in de Sitter spacetime and discuss how the classical nature of quantum fluctuations 
appears. Section V is devoted to summary. We use units in which c = h = G = 1 throughout 
the paper. 



II. TWO DETECTORS SYSTEM 



We consider a system with two equivalent detectors interacting with the massless scalar 
field in an expanding universe [ll|, 13|. The detectors have two energy level states | t)? I i) 
and their energy difference is given by Q. The interaction Hamiltonian is assumed to be 



V = g{t){a\ + a-^)(l){xA{t)) + g{t){a+ + (y-s)(t>{xB{t)) 
where ct+, are raising and lowering operators for the detector's state: 

|t)ai, a- = |i)(t|. 



(1) 



Two detectors are placed at XA^sif) and {t,XA,Bif)) represent their world lines. We assume 
that the detectors are comoving with respect to cosmic expansion. Strength of the coupling 
is controlled in accord with the following Gaussian window function 



g{t) = 5-0 exp 



{t - to)- 
2a2 



(2) 



This window function approximates the detector being "on" when |t — to| ~ o" and "off" the 
rest of time. We assume that the detectors are both down state initially {t — )■ — oo) and the 
scalar field is vacuum state |0). Thus, the initial state of the total system is |\l/o) = | ii)|0). 
Then, in the interaction representation, the final state {t — ?■ +oo) of the total system after 
interaction becomes 



dtiVi - - 
2 



dt,dt2T[ViV2\ + 



l^o) 



(^1 - It [$;^$+ + <i.^<i.+]^ I ii)|o) -^$+1 ti)|o) - ^$+1 it)|o) 



(3) 



-T + $+$+] I tt)|0) + O(^?^ 



where T is time ordering and symbols with a tilde denote quantities in the interaction 
representation. We defined a field operator 



dt,g{ti)e^'''''4>{t,,XA,B{ti))- 
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As we are interested in the detectors' state, by tracing out the degrees of freedom of the 
scalar field, the state for the detectors' degrees of freedom becomes 



P 



/X4 X \ 

E Eab 
Eab E 
\X* 1-2E-Xij 



(4) 



where we use the basis {ft? ti? it? ii] for this matrix representation of the state. The matrix 
elements of the state (HI) are 



/oo nt\ 
dti / dt2gig2e'''^''^''\(t>{ti,XA)(t>{t2,XB)), 
-00 J —00 



E 



AB 



X. 



'OO 

00 



dti / dt2gig2e 



(0(ti, a;A)0(t2, Xb)), E = EAB{r = 0), 



dta 



-00 



dt^gsgi 



X 



dti 



C^t2^l(72e-^^(*«+**)e^^(*i+*^)(0(t4, XB)<P{t3, XA)<P{tl, XA)<P{t2, Xb)). 



In X4, the four point function can be written using two point functions 

(0(^4, XB)(p{t3, XA)(t>(ti, XA)(p{t2, Xb)) = {<P{ti, 03^)0(^1, OJa)) (0(^3, XA)<P{t2, Xb)) 

+ (0(^4, XB)(j){t3, XA)){(f){tl, XA)(j){t2, Xb)) 
+ {(p{t3, XA)(pitl, XA)){(p{ti, XB)(t){t2, Xb)). 

E,Eab are 0{g^) and X4 is O^g"^) quantities. For the purpose of obtaining entanglement 
only, it is not necessary to evaluate X4. We need X4 to obtain quantum mutual information 
in which quantum part of correlations are encoded. By changing integration variables to 



tl +to 



t1 -t. 



+ tA 



t-i — td 



Xi 



yi 



X2 



1/2 
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we have the following expressions for the matrix elements: 

/oo POO 
-OO Jo 
poo poo 

Eab = 2gl / dxe-^"'"' / dye-y'/^'-^'^'y D+{x + to, 2/, r), 



oo 
oo 



/oo /"OO /"OO /"OO 

-oo Jo ./-oo Jo 



(5) 



X 



n+ , a;i + 1/1 + X2 - ^2 , , X2-y2-xi- yi 
^ I r to, 7, 1^ 



X D- 



xi-yi + X2 + y2 , X2 + y2-xi + yi 



+ t 



0, 



D^{xi + to, yi, r)D+{x2 + to, 2/2, 
, n+ / +^1 + ^2 + 2/2 , , X2 + y2~xi-yi 

+ ^ I ~ r to, ,U 



2 
X D 



+ I xi~yi + X2-y2 , , X2-y2- xi+yi 
I 2 + ^0' 2 ' 







where we have introduced the Wightman function for the scalar field 

D^{xi,yi,r) = {(j){ti,XA)(i){t2,XB)), r =\xa- Xb\. 

III. ENTANGLEMENT AND CORRELATIONS OF DETECTORS 

As the response of detectors due to interaction with the scalar field is given by the state 
(jl]), we can extract information on entanglement and correlations of the scalar field indirectly 
by analyzing this state. 



A. Entanglement of detectors 

As a measure of the entanglement between two detectors, we consider the negativity 



14l | defined via a partially transposed operation to the density matrix with respect to 



detector B's degrees of freedom. The eigenvalues of the partially transposed density matrix 



are 



A = E±\X\+0{g^). 

The negativity is defined using the eigenvalues of the partially transposed density matrix 



max 



OAX\-E 



Ai<0 



From here on, we designate the following quantity as the negativity 

Af= \X\-E. 



(6) 



(7) 



The negativity gives the necessity and the sufficient condition of the entanglement for two- 
qubit systems jl, 0]. Thus two detectors are entangled when J\f > and separable when 
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A/" < 0. For separable initial states of detectors A/" < 0, A/" > after interaction with the 
scalar field implies the scalar field is entangled because entanglement cannot be generated 
by local operations. 



B. Correlations of detectors 

Using Bloch representation, the state (jlj) can be written as follows 
p^^ = ^\l'S)I + a-(T®I + I®b-cr + q^ct^ ® cr„ I , (8) 



1 f ^ 
-ll^I + a- (T®I + I®b-cr + ^ ( 

\ e,m=i 



/2{Eab + Xr) -2Xj 
a = b= {0,0, -1 + 2E + 2X4), = -2X, 2{Eab ~ Xr) 

\ 

where / is the identity operator, {cri, cr2, (T3} are the Pauli spin matrices, Xr = Re(X) and 
Xj = Im(X). To quantify quantumness of quantum fluctuations in de Sitter spacetime, 
we want to consider the "classical" and "quantum" part of the correlation between two 
detectors. The "classical" part of correlation is defined through a local measurement on each 
detector. By a measurement here we mean the von Neuman type; complete measurement 
consisting of orthogonal one- dimensional projectors. 



1. Classical mutual information 



To define the classical part of the correlation between two detectors (two qubits), we 
perform a local projective measurement of detector states. Of course, it is not possible 
to perform measurement of the scalar field in the inflationary era directly. We consider 
the following measurement procedure as a gedanken experiment to explore the nature of 
quantum fluctuation. The measurement operators for each detector are 



M. 



± 



M 



2 > --^ 2 ' |nA| = |nB| = l (9) 

where d= denotes output of the measurement, n^, ris represent the internal direction of 
measurement. The joint probability pjk attaining measurement result j for detector A and 
k for detector B (j, k = ±1) is obtained as 

tB AB\ 



p,k = tT{Mf p 
1 



4 



(1 + (j>.n^)(l + {k)hnf) + ijk) Y,{cim - 

Im 

2(-Eab + Xji) —2Xi 
— 2X/ 2{Eab — Xr) 



ebm)nfn^ 








(10) 











-4^^ + 4x4^ 



The probability pj attaining a result j for detector A and pk attaining a result k for detector 
B are 

k j 
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Using the joint probability (fTOj) obtained by the measurement, the classical mutual infor- 
mation Ic is defined by [l5| 

Ic{p) = H{pj) + H{pk)- H{pjk), H{pj) = pj log2 Pj, H{pjk) = - Pjk log2 Pjk 

j jk 

where H{p) is a Shannon entropy for a probability distribution p. The explicit form of Ic{p) 
using Pjk is 

/c(p)^Ef.^'°^.te)- (11) 

jk 

The measure of classical correlation is defined via the maximization done over all possible 
projective measurements:^ 

C{p)= sup Icip). (12) 

{nA,nB} 

We evaluate the classical mutual information of pjk separately for the cases where the 
vectors Ua^b are parallel to the z axes or not. For the case where ua and ub are both 
parallel to the z axes, we have 

Icip) = ^{e'-X. + X, In ) + Oig% (13) 

For the case where one of ua or is parallel to the z axes, 

Icip) = 0{g'). (14) 
For the case where both ua and are not parallel to the z axes. 



1 \T.em=x,yCim{nA)e{nB) 



2 

,6\ 



_ '^{Eab + Xr) —2Xi 



By introducing new direction vectors 



"■A = ^} , I ^^"^ I , riB 



nAx \ - ^ nsx 



+ {riAy)' V^yJ v^(^^^^FTK^ y'^y' 

we obtain ^ 

Ic{p) = I XI cem{nA)e{nB)m J • (16) 

\ / 

As the eigenvalues of the matrix c is 2{Eab ± \X\), the maximum of classical mutual infor- 
mation for this case is given by 

Ic,mUP) = :^{\EAB\ + \X\y. (17) 



^ This definition of tlie classical correlation is based on two-sided measurements of correlations. We can 
also define the classical correlation using one-sided measurements of correlations. 
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Thus, the classical correlation of the state (|8]) obtained by the local projective measurement 
is 



max 



1 



In 2 



X4 + X4ln( ^ ) ), -^[\Eab\ + \X 



E^ 



In 2 



(18) 



We comment on the relation between the classical mutual information and a correlation 
function of qubit variables. By the local projective measurement of the state (|8]), we obtain 
the following expectation values for qubit variables: 

{ua ■ (t) = a - nA, {riB ■ cr) = b ■ ub, 
ijiA ■ a ®nB ■ cr) =^ Cem{nA)e{nB)m- 



By introducing a fluctuation part of qubit variables as Acr^ 
function for fluctuations of qubit variables is 



(Aan^Adng) = ^(Qm - athm){nA)t{nB)m = ^Qm(^A)K^i?)m 

£,m 



n - (T — (n cr), the correlation 

(19) 



where we have assumed n^, are not parallel to z axis at the last equality. Therefore, the 
classical mutual information flTBl) for joint probability pjk corresponds to the square of the 
correlation function of the fluctuation part of the qubit variables. 



2. Quantum mutual information and quantum discord 

The quantum mutual information Iq of the bipartite state p^^ is defined independently 
of measurement procedure: 



5(P") + Sip"") - Sip 



^AB\ 



S{p) = -tr(plog2p) 



(20) 



where S{p) is the von Newmann entropy for the state p. This quantity represents the 
total correlations of the bipartite system including both the quantum and classical parts of 
correlations. 

We evaluate the quantum mutual information for the state ([8]). The eigenvalues of the 
state ([8]) are 1 — 2E + |Xp — X4, E ± Eab, — The reduced density matrix for the 

subsystem A is 



E + X^ 
1-E-Xi 



Thus, the quantum mutual information for the two detector system is given by 



Up) 



1 

h72 

+ 



-2E In E+{E + Eab) HE + Eab) + {E - Eab) HE - Eab) 



(21) 



ln2 L 



^2 -X4 + -2X4lnE + (X4- |Xnin(X4 - |X| 



The quantum discord 1^ 18| is introduced as the difference between the quantum mutual 
information and the classical mutual information 



Q{p) = Hp)-c{p). 



(22) 
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For arbitrary local projective measurements, it can be shown that Q > 0. Thus, the quantum 
mutual information can be decomposed into a positive classical mutual information and a 
positive quantum discord. The necessity and sufficient condition for zero quantum discord 
is expressed as the state has the following form 1^ : 



P^^ = Y.P^kMfM^. (23) 

For this state, measurement ([9]) does not alter the form of state because M^^, are projec- 
tion operators. In this sense, the state with zero quantum discord can be termed a classical 
state and the quantum discord represents the quantum part of the total correlations. For 
pure state, we have Q = C and the state with zero quantum discord has no classical corre- 
lations. 



3. Bell- CHSH inequality 

Related to classical correlations obtainable via local measurements, we consider the ques- 
tion whether correlations derived under the state ([S]) admit a local hidden-variable (LHV) 
description; measured correlation functions can be mimicked by classical distribution func- 
tions. Let us consider the following operator (Bell operator): 

ScHSH = a-(T®{b + b')-(T + a' ■(T^{b-b')-(T (24) 

where a, a', b, b' are real unit vectors. Then, the Bell-Clauser-Horne-Shimony-Holt(CHSH) 
inequality |20| is 

|(Schsh)|<2. (25) 

If the state admits a LHV description of correlations, then this inequality holds. Violation of 
this inequality means existence of nonlocality. The two qubit state violates the Bell-CHSH 
inequality if and only if the following condition is satisfied [2l| : 

M(p) > 1, M(p) = the sum of the two largest eigenvalues of the matrix c^c (26) 

where the matrix Cjk = ti^paj cTfc). For the state ([8]), 

(Cll Ci2 \ 
ci2 C22 , cn = 2{Eab + Xr), 
C33/ 

C22 = 2{Eab - Xr), C12 = -2Xj, C33 = 1 - 4E. 

As eigenvalues of Cjk are 1 — 4:E,2{Eab ± |X|), the sum of the square of the two largest 
eigenvalues of c cannot exceed unity and the Bell-CHSH inequality holds. However, this does 
not mean the LHV description of correlations is possible; holding the Bell-CHSH inequality 
is only a necessary condition for the LHV description and does not guarantee existence of a 
LHV. Actually, by passing each detector through the filter 

^-^■^ = (0 !!) ' 
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there is a possibility revealing hidden nonlocality of the state |22|. After passing through 
the filter, the state is transformed according to p — )■ p' = (/^ ® /_b)p(/a ® fs)- Assuming 
that = 0{g'^)^ the transformed normalized state is 



P' 



( X4 




trp' X4 + 2^^E + rf 
and the matrix c is transformed to 



2x* 



X4 + 27^2^ + r/4 









rfX\ 






















1 


-Xi 








/ -£"715 + Xr 

Xi Eab — Xr 

X4 - 27^2^ + r/^ 



v 











27^2 



Assuming that the state is entangled = |X| — > 0, eigenvalues of c' are 

2rf\X\ AErj'^ 

^4 + 774' ~ xTT^ 
and the condition for violation of the Bell-CHSH inequality M{p') > 1 is 

v' - ^-§^v' + X, < 0. 



For existence of a real rj satisfying this inequality, we need [10 

IXI^ > I6X4EI 



(27) 



(28) 



(29) 



This provides a sufficient condition for violation of the Bell-CHSH inequality and existence 
of a hidden nonlocality. 
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IV. BEHAVIOR OF ENTANGLEMENT AND CORRELATIONS OF SCALAR 
FIELD 



The matrix elements Eab , X, X4 of ([5]) can be evaluated using numerical integration 
after considering contributions of poles in integrands by contour integration on a complex 
plane 23| . However, numerical integrations of these functions are not so easy especially eval- 
uating the fourfold integral of X4. Hence, in this paper, we consider asymptotic estimation 
of these functions to derive analytic approximate forms. We assume parameters Qa ^ 1 
with other dimensionless combinations of parameters such as HQa^ contained in are 
kept order unity. That is, we are considering the asymptotic behavior of Eab, X, in the 
range of parameters 

fia > 1, Ha 

where H is the Hubble parameter. By rescaling the integration variables x and y, 

/OO roo 
-00 Jo 

/OO POO 
-00 J —00 

For Qa ^ 1, x and y integrals in X can be evaluated approximately at the saddle point 
X = i and ?/ = of the integrand. After performing Gaussian integrals about these saddle 
points, we have 

/OO roo 
-00 Jo 

In the same way, expression for Eab is 

/OO POO 
'OO J —00 

Thus, the asymptotic forms of X, Eab for Qa ^ 1 are obtained as 
X ^ -27Tgla^e-^^''^'D+{to + i^a^ 0, r), 

Eab ~ 27r(7oVe-(^'^)'Z}+(to, -^^^a^r), (30) 



E' + El^ + \X\'. 



9o^ 

~ ^AB 

The negativity is 

M ^ 27r^^^2g-(f^a)2 \\jj+^t^ + iQ^2^ 0, r) I - D+{to, -i^a^, 0)1 . (31) 

Using the asymptotic form of X4, we can reduce the condition of hidden nonlocality f l29|l . For 
the entangled region |X| ^ E, using the asymptotic form of X4 ^ E"^ + |Xp + E]^^ I^P) 
the condition of hidden nonlocality f l29|) reduces to 

|X| > AE (32) 
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and provides a stronger condition than that for entanglement |X| > E. 

Now, let us we consider entanglement and correlations of the scalar fields in Minkowski 
and de Sitter spacetime. 



Minkowski vacuum 



The Wightman function of the massless scalar with Minkowski vacuum state is 



47r2 A{y - it) 



(33) 



where a positive small constant e > is introduced to regularize ultraviolet divergence in 
the Wightman function (see Appendix). Using f l30|) . 



X = -{2-1X91) 



1^ 



E 



AB 



47r2 (r/a)2 + 4(fi(T)2" 



The negativity is 



Af 



and the = line is given by 



9o_ -n^a^ 
8n 



a 



(na) 



2a' 



(34) 



(35) 



For r/a < 2Qa, the negativity is positive and detectors are entangled (Fig. [T]). As the 
prepared initial state of detectors is separable and entanglement cannot be created via 
local operations, this entanglement is due to the scalar field and it is swapped to a pair 
of detectors through interaction. For large separation r/a > 2Qa, detectors are separable. 
However, this does not mean the scalar field is separable for this separation; by increasing 
detector's parameter a with fixed r, it is possible to make A/" > 0. For any separation r, we 
can always choose an appropriate value of a that makes detectors entangled. Thus, we can 
say that the Minkowski vacuum is always entangled. 
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0.0 



0.5 



1.0 



2.0 



FIG. 1: Contours of negativity for the Minkowski vacuum. Two detectors are entangled with 
parameters in the shaded region. Darker areas correspond to larger negativities. In the region 
above the dotted line, the Bell-CHSH inequality is violated after operation of local filtering. 

Now, let us consider behavior of correlations. The left panel in Fig. |2]shows r dependence 
of EjHb- They are monotonically decreasing functions with respect to r. In the entangled 
region r < 2Qa'^, Eab is nearly constant E^b ~ E and in the separable region r > 2ficr^, 
Eab ^ \X\ < E. 



0.1 




rn 



FIG. 2: r dependence of correlations (crO = l,gQ = 0.1). The right panel shows r dependence of 
Iq and C. In both panels, dashed lines represent negativity. 

r dependence of the quantum mutual information Iq and the classical correlation C is shown 
in the right panel in Fig. [21 These correlations also decrease monotonically with respect to r. 
In the entangled regions, the difference of these two correlations (quantum discord) increases 
as r increases until two detectors become separable. In the separable region, quantum discord 
remains constant. 
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50.0 



FIG. 3: r dependence of the ratio Cjlq (crO = 1, 50 = 0.1). The dashed hne represents negativity. 



As shown in Fig. |3], the ratio of classical correlation to total correlation Cjlq is a good 
indicator of separability for this system. In the entangled region, the ratio decreases mono- 
tonically and in the separable region, the ratio approaches a constant value. We can derive 
this behavior using the following r dependence of X, Eab'- 



r <^ Tc : 
r ^ Tc : 

where Tc = 2Qa'^. For r <^ r^, 

C^max<! — iXPln^, -^IXI 
In 2 ' ' E ' In 2 ' ' 



Eab ~ -E" -C |X| (entangled region), 
\X\ ^ Eab ^ E (separable region). 



(36) 



In 2' 



ixrinE. 



(37) 



In this region, as |X| ^ E, the ratio is 



ln|X| 
InE 



and this reproduces decreasing behavior of the ratio with respect to r. For r ^ Tc, 

C - ^El,, Iq ^ (38) 
and the ratio C/Iq approaches a constant value independent of r and its value is given by 



C 

— ^AE 



2' 



(39) 



The behavior of r dependence of the ratio C/Iq changes at |X| E and this corresponds 
to the separable line M = The separability condition A/" < comes from the positive 
partially transposed criterion of the state [IJ]. Although the quantity C/Iq is introduced 
independently of the positive partially transpoed criterion, the change of its r dependence 
ai M ^ Q arises and it corresponds to the separability condition. 
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B. Scalar field in de Sitter spacetime 

In de Sitter spacetime with a spatially flat slice 

the Wightman function of a massless conformal scalar field with the conformal vacuum 
state l?! is 

DLf = [- sinh^(iJ(y - ^e)) + e''^%Hr/2y] . (40) 

The Wigthman function of a massless minimal scalar field with the Bunch-Davis vacuum 
state is ^ 

^iin = ^inf + ^2, (41) 
tt2 

D2 = [Ei [i{Hr - 2e-^^sinh(/7(y - ie)))] + Ei [i{-Hr - 26"^^ sinh(i/(|/ - te)))]] . 

The massless minimal scalar field in de Sitter spacetime suffers from infrared divergence. 
Hence we have introduced the infrared cutoff = H m D2 which corresponds to a maximal 
comoving size of the inflating region; r is the comoving distance between two spatially 
separated points and satisfies r < H~^. We have assumed that inflation starts at t = and 
the comoving distance r must be smaller than the size of the inflationary universe H^^ at 
t = 0. The asymptotic form of X, Eab are 

-^conf = EjiBconl = -^ABl, (42) 

^min = Xi + X2, Eab min = EaBI + EaB2 



where 



X, = -(2.gl) \^^ )Z\ . 



0)' 

Eabi = i2ngl) \^^ ^ . , , (43) 



47r2 {Hrpj 



2 

47r2 Asin^{Hna^) + {Hrp^ 



X2 = {27iglf-^{Haf {Ei (-ze'^^^ifr,) + Ei (^e'^^^ifr,) } , 

Eab2 = -{2T:gl)^-^{Haf\Ei {-ie-"''Hrp + 2e-^*« siniSlHa^)) 

+ Ei {ie~"^°Hrp + 2e-"^° sm{nHa^)) | 

where Hto is the e-folding time at which detectors are switched on. The physical distance 
between detectors at this instance is Vp = e^^°ro with tq < H~^. Unlike the Minkowski vac- 
uum case, physical distance between two detectors increases in accord with cosmic expansion 
and entanglement between detectors changes in time. 
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FIG. 4: Contours of negativity for the conformal invariant scalar field (left panel) and the minimal 
scalar field (right panel) (0,/H = l,ro = 0.1). Two detectors are entangled with parameters in 
the shaded region. In the region above the dotted line, the Bell-CHSH inequality is violated after 
operation of local filtering. 



For sufficiently large e-foldings Hto ^ 1, 



-(27r(7, 



E 



AB rain 



E 



AB2 



47r2 

g-(0(7) 



{Hcxy[-ln{Hrp) + Hto], 



47r2 



Asm' (Hna^) 



and the negativity for the minimal scalar field is 



2n 

The separability condition A/" < yields 

Hrp > 2sm''{HQa')exp 



Hr-n 



1 



Asm' HVta' 



Asm\Hna') 



1.0, 



(44) 



(45) 



where the numerical value is obtained for HQa^ = 1. Comparing this to the Minkwski 
vacuum case, we can find sufficiently large at which detectors are separable for any value 
of a detector's parameters fl, a (see Fig. 4). As Vp grows in time to, two entangled detectors 
Tp < evolve to a separable state after their separation exceeds the Hubble horizon 
scale. This behavior is consistent with our previous analysis of entanglement using the 
lattice model and the coarse-grained model of the scalar field jH, 0]: bipartite entanglement 
of the scalar field between spatially separated regions in de Sitter spacetime disappears 
beyond the Hubble horizon scale. For the conformal invariant scalar field, the separability 
condition becomes 



Hr„ > 2 sin HVla^ ~ 1.4. 



(46) 
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The numerical value is for HQcr^ = 1. Contrary to the minimal scalar field, the entangled 
region extends to the super horizon scale. Although the size of the entangled region is larger 
than that of the minimal scalar field, the two detectors are separable when their separation 
reaches scale (1461) and this behavior is the same as the minimal scalar field. 

For the entangled state, we can check the violation of the Bell-CHSH inequality. The 
condition of violating the Bell-CHSH inequality after operation of local filtering is 



Hrp < sm{Hna^) ~ 0.84 



1 



sin^ Hna"^ 



0.34 (minimal scalar), 
(conformal scalar). 



(47) 



Thus, for both scalar fields, the appearance of hidden nonlocality is possible only on the 
subhorizon scale and violation of the Bell-CHSH inequahty cannot be detected on the super 
horizon scale. 

Behavior of correlations is completely different for these two types of scalar fields. Figure E] 
shows Tp dependence of functions for the conformal invariant scalar field and the 

minimal scalar field. 
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FIG. 5: rp dependence of \X\ and Eab (left panel: conformal invariant scalar field, riglit panel: 
minimal scalar field, Ha = 0.5, ^l/H = l,ro = 0.1, go = 0.1). rp is related to to as rp = roe^*°. 
The dashed line represents the negativity. 



For the conformal invariant scalar field, these functions decay as ~ on super horizon 
scale, which is the same behavior as the Minkowski vacuum case. On the other hand, for 
the minimal scalar field, these functions approach a constant value on the super horizon 
scale due to accumulation of long wavelength modes of quantum fluctuations. The different 
behavior of X, Eab leads to different behavior of classical and quantum parts of correlations 
of these scalar fields. 
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In Fig. ini the classical and total correlation C, Iq of the minimal scalar field also approach 
constant values on the super horizon scale while these correlations of the conformal invariant 
scalar field decay as ~ r~'^. 




0.5 1.0 5.0 10.0 50.0 0.5 1.0 5.0 10.0 50.0 



FIG. 7: Tp dependence of the ratio C/Iq (left panel: conformal invariant scalar field, right panel: 
minimal scalar field, Ha = 0.5,Q/H = l,ro = 0.1, ^fp = 0.1). The dashed line represents the 
negativity. 



In Fig. [TJ the ratio C/Iq for the conformal invariant scalar field approaches constant for 
rpH ^ 1 and this behavior is the same as the Minkowski vacuum case (see Fig. 3). However, 
for the minimal scalar field, value of this ratio has Inr^ dependence and increases as Vp 
increases: 

For a sufficiently large distance (large e-folding Hto) given by 



roj gl{Hay ~ gl \H 

C ^ Iq and this means the quantum state of detectors approaches the zero quantum dis- 
cord state and can be regarded as classical in the sense that measurement procedure does 
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not alter the state. As discussed in Ref. 12|, classical and quantum channel capacities of 



communication via detectors beyond the light cone identically vanishes and we presumably 
expect these capacities to be zero beyond the Hubble horizon. Thus the correlation with 
zero discord on the super horizon scale originates from quantum fluctuations of the scalar 
field. Therefore, this behavior of correlations supports the long wavelength quantum fluctu- 
ations of the massless minimal scalar field in de Sitter spacetime being treated as classical 
fluctuations. 



V. SUMMARY 

We investigated quantum and classical correlations of the quantum field in de Sitter 
spacetime using the detector model. Entanglement of the scalar field is swapped to that 
of two detectors interacting with the scalar field and we can measure the entanglement of 
the quantum field by this setup of experiment. In de Sitter spacetime, the entanglement 
between detectors disappears on the super horizon scale and this behavior is consistent with 
our previous analysis using the lattice model and the coarse-grained model of the scalar 
field [2I . However, the behavior of correlations shows different behavior depending on the 
type of scalar fields. For the massless minimal scalar field, the ratio of classical correlation to 
the total correlation approaches unity for sufficiently large e-foldings. On the other hand, for 
the massless conformal scalar field, that ratio approaches a constant value smaller than unity 
and the condition for classicality is not achieved. These results support the long wavelength 
quantum fluctuation of the minimal scalar field being treated as classical fluctuations and 
becoming seed fluctuations for the structure in the our Universe. 

As an application of our analysis presented in this paper, it is interesting to consider 
quantum effects in analogue curved spacetimes proposed using Bose-Einstein condensates 
or ion traps j^^. In these experiential setups of analogue models, we can directly measure 
entanglement and classical and quantum correlations of quantum fluctuation using detectors 
in the laboratory. We expect that investigation in this direction will increase understanding 
of the quantum and classical nature generated during the inflation. 
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Appendix A: Wightman function of massless minimal scalar field 

In de Sitter spacetime with a spatially flat slice, the massless minimal scalar field obeys 
the following equation of motion 

+ 3H^ - e"^^* VV = 0. (Al) 
The quantized field with the Bunch-Davis vacuum state is 



x) = j ^^-^ [m d, + r,it) a_,t) e^'^-, = ^ _ 1 ) e--'' (A2) 
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where rj — —e H is conformal time. The Wightman function is 



1 r°° 

^ dkk'jo{kr)Mr),)f^{r]2), 



27r2 



e-'^'smkre-'''^" + 



I 



Xi - X2\ 

dke'"^ sin kr 



d_ 
dk 



-ikArj 



(A3) 



;A77 = 77i -772, 



where we have introduced a damping factor e with a small positive number e to regularize 
ultraviolet divergence of the k integral. The first integral is 



y 



(A4) 



This term is the same as the Wightman function for the massless conformal invariant scalar 
field. The second integral diverges at A; = 0, hence we introduce a lower bound ko of k 
integral: 

r 

-Ei {—ko{ir — iArj — e)) — Ei {—ko{—ir — iArj — e)) 



Do 



(A5) 



As the value of /cq, we choose k^ = H which corresponds to the size of the universe (horizon 
scale) at the beginning of inflation. 
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